A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
08(09.04.1;09.07.1;13.21.3) 



ASTRONOMY 
AND 
ASTROPHYSICS 
5.2.2008 



G\ 
ON 

> 
o 

in 



> 
o 

(N 



ON ■ 

Oh; 

6 ■ 

a: 



Constraints on the properties of the 2175A interstellar 
feature carrier 

F. Rouleau, Th. Henning and R. Stognienko 

Max-Planck-Gesellschaft, AG "Staub in Sternentstehungsgebieten" , Schillergafichen 2-3, D-07745 Jena, Germany 
Received 2 July 1996 / Accepted 15 November 1996 



Abstract. Constraints on the possible shape and clus- 
tering, as well as optical properties, of grains responsible 
for the 2175 A interstellar extinction feature (interstel- 
lar UV bump) are discussed. These constraints are based 
on the observation that the peak position of the inter- 
stellar UV feature is very stable (variations 1%), that 
the large variations in width (^ 25%) are uncorrelated 
with the peak position except for the widest bumps, and 
that the shape of the feature is described extremely well 
by a Drude profile. The UV extinction of small graphite 
grains is computed for various clustering models involv- 
ing Rayleigh spheres. It is shown that compact clusters 
qualitatively satisfy the above observational constraints, 
except that the peak position falls at the wrong wave- 
length. As an alternative to graphite to model the optical 
properties of the interstellar UV feature carrier, a single- 
Lorentz oscillator model is considered, in conjunction with 
a clustering model based on clusters of spheres. Intrinsic 
changes in the peak position and width are attributed to 
variations in chemical composition of the grains, impact- 
ing upon the parameters of the Lorentz oscillator. Fur- 
ther broadening is attributed to clustering. These models 
are shown to satisfy the above observational constraints. 
Furthermore, the correlated shift of peak position with 
increased width, observed for the widest interstellar UV 
features, is reproduced. Models involving a second Lorentz 
oscillator to reproduce the FUV rise are also considered. 
The impact of this extra Lorentz oscillator on the peak 
position, width, and shape of the bump is investigated. 
Synthetic extinction curves are generated to model actual 
ones exhibiting a wide range of FUV curvatures. Physical 
mechanisms which might be of relevance to explain the 
variations of these optical properties are discussed. 
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1. Introduction 

Ever since its discovery, the prominent 2175 A UV feature 
in the interstellar extinction curve has been attributed to 
graphite - or at least some form of carbon-rich material 
- due to cosmic abundance constraints, and the fact that 
small graphite particles produce a feature in roughly the 
right wavelength range. 

The interstellar extinction curves along different lines 
of sight can be reproduced amazingly well by the following 
simple parameterization (Fitzpatrick & Massa 1986; 1988; 
1990) 
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where x is the inverse wavelength (wavenumber in /mi -1 ), 
A is the magnitude of extinction, a; max is close to the peak 
wavenumber of the UV feature, and ai, a-2, 03, 04 are 
fitting parameters. The far-UV curvature is described by 
/ PUV (x) = 0.5392(x - 5.9) 2 + 0.0564(a; - 5.9) 3 for 5.9 < 
x < 8.0 /mi -1 and fFw(x) = for x < 5.9 /trn -1 . The 
function \{x — x^ ax /x) 2 + 7 2 ] -1 is referred to as a "Drude" 
profile, though it is a profile produced by small spheres 
described either by a single Lorentz or Drude oscillator 
(Bohren & Huffman 1983). This fit reproduces most if not 
all interstellar extinction curves known to date, with no 
systematic deviations or additional features or shoulders 
(at least, none beyond the observational errors). 

The main observational constraints concerning the in- 
terstellar UV feature (term with coefficient 03 in Eq. [1]) 
are (Fitzpatrick & Massa 1986; Jenniskens & Greenberg 
1993): 

1. the remarkable constancy of its peak position, ir max = 
4.60 ± 0.04 /mi -1 , though its small (Si 1%) variations 
are larger than observational errors. 

2. the wide range of variations in its width, 7 = 1.0 ± 
0.25 /mi -1 (i.e. <> 25%). 
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3. the fact that the variations in peak position and width 
are uncorrelated, except for the widest bumps, i.e. 7 k5 
1.2 fixa. , for which a systematic shift to larger x is 
observed (Cardelli & Savage 1988; Cardelli & Clayton 
1991). The lines of sight for which this is observed all 
pass through dark, dense regions. 

Various lines of sight show peculiar extinction curves 
depending on the type of environment they pass through. 
For example, in the hydrogen poor circumstellar envi- 
ronment of R CrB, the bump is weaker and shifted to 
4.17 /jitT - 1 (2400 A). Around carbon-rich (and hydrogen- 
rich) asymptotic giant branch stars, the bump is consid- 
erably weakened or absent (e.g., Snow et al. 1987). In Hn 
regions, the relative strength of the bump (parameter 03) 
is weaker than in most other environments (Jenniskens & 
Greenberg 1993). However, the range of bump widths is 
similar to that of the diffuse interstellar medium. Broader 
bumps in Hn regions have their peak position shifted to 
smaller x. In other dense environments, like "bubbles" 
(regions showing loops and filaments characteristic of ma- 
terial swept up by stellar winds of OB stars in regions of 
recent massive star formation), the strength of the bump 
is similar to that of the diffuse ISM and no correlation is 
observed between its width and its peak position. 

In this paper constraints on the optical properties of 
the purported interstellar UV feature carrier are derived 
considering rather general arguments related to chemical 
composition (as modelled by graphite or Lorentz oscillator 
models) and clustering (based on direct computations and 
interpretation using a spectral representation). In Sect. [2] 
some current theoretical models attempting to explain the 
characteristics of the interstellar UV feature are discussed. 
Their main shortcomings are emphasized. In Sect, [j] the 
profile of the UV feature of graphite is computed for var- 
ious arrangements of touching spheres. The results are 
compared qualitatively to the observational constraints. 
Apart from the peak position falling at the wrong wave- 
length, the qualitative agreement is found to be good for 
compact clusters. But variations in chemical composition 
cannot be included directly in such a model. Therefore, in 
Sect. H a series of single-Lorentz oscillator models are used 
in conjunction with clustering to investigate the range of 
parameters (which simulate variations in chemical compo- 
sition) consistent with the observational constraints. The 
clustering is modelled via a spectral representation formal- 
ism. In Sect. |^ the effects of adding a second Lorentz os- 
cillator to explain the FUV rise are discussed. Interstellar 
curves along specific lines of sight exhibiting a wide range 
of FUV curvatures arc modelled and plausible (though 
non-unique) optical constants for the bump grains along 
these lines of sight are derived. In Sect. ^| physical mech- 
anisms relevant to the expected optical properties of the 
UV feature carrier are discussed, in particular, dehydro- 
genation and UV processing. 



2. Current theoretical models 

Detailed theoretical comparisons between electromagnetic 
scattering models involving small graphite grains and the 
interstellar UV feature are complicated by the fact that 
graphite is a highly anisotropic material. In particular, an 
approximation must be used when computations involv- 
ing Mie theory (valid only for isotropic spheres) are carried 
out. The MRN model (Mathis et al. 1977) and its variants, 
using a size distribution of spherical grains of silicate and 
graphite, can explain the mean extinction curve, but fail to 
satisfy the above observational constraints. Draine (1988) 
has studied the UV feature produced by graphite parti- 
cles using the discrete dipole approximation (DDA). This 
method is ideally suited for handling the anisotropic di- 
electric tensor of graphite. He found that only particles of 
small elongation with equivalent radii in the range 100- 
200 A could provide a reasonable fit to typical interstellar 
UV features. 

Draine & Malhotra (1993) studied the variations in 
peak position, width, and strength of the bump for various 
models based on a size distribution of graphite grains to 
see whether they were compatible with the observational 
constraints. DDA calculations included spherical graphite 
grains with an ice coating, spheroidal graphite grains, and 
graphite spheres in contact with silicate spheres. All mod- 
els using variations in shape or coating produced corre- 
lations between the peak position and width. Therefore, 
they concluded that the variations observed must be due 
to changes in the dielectric properties of the grains either 
through impurities or surface effects, rather than purely 
"geometric" effects. 

Mathis (1994) considered a model consisting of a 
graphite oblate spheroidal core and a coating of mate- 
rial represented by an appropriately chosen single Lorentz 
oscillator. The weakness of this model lies in the fact that 
the shape of the grain had to be unreasonably fine-tuned 
in order to reproduce the stability in peak position of the 
UV feature. Having a coating to broaden the bump and 
eliminate correlations between its width and its peak po- 
sition can be considered a second order effect, the overall 
shape of the grain being the first order effect. Therefore, 
implicit in this model is the unlikely assumption that all 
interstellar grains along every possible line of sight have 
exactly the same shape. Any deviation in shape shifts the 
peak position outside the observed range. Furthermore, 
variations in shape introduce a correlation between the 
strength, the width, and the peak position of the feature. 
Since the width and peak position of the narrowest in- 
terstellar UV features are narrower and shifted to smaller 
wavenumbers, respectively, compared to those of Rayleigh 
graphite grains, Mathis had to "tinker" with the dielec- 
tric function e± of graphite tabulated by Draine (1985). 
This is somewhat justified in view of the lack of agreement 
between laboratory measurements carried out by various 
authors in that spectral region (see e.g., Draine & Lee 
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1984), but it introduces additional free parameters. Note 
that virtually all of these graphite optical constant mea- 
surements yield a plasmon resonance for Rayleigh spheres 
at wavenumbers significantly larger than 4.60 /Ltm -1 . 

Henrard et al. (1993) considered a model in which 
the bump carrier was assumed to consist of small spheri- 
cal onion shells of graphite. Again, an unreasonable fine- 
tuning in the number of shells was required to reproduce 
the peak position of the interstellar UV feature. 

3. Clustering effects on the graphite UV feature 

In this section we investigate whether clustering of grains 
composed of graphite can lead to an increase of the width 
of the UV feature without appreciably changing its peak 
position. 

Draine (1988) and Draine & Malhotra (1993) have con- 
firmed the good agreement between the so called "1/3- 
2/3" approximation and DDA computations taking ex- 
plicitly the anisotropy of graphite into account. In this 
approximation extinction cross sections C(e\\) and C(e±) 
are obtained separately for isotropic particles having the 
components of the dielectric tensor parallel and perpendic- 
ular to the c-axis of graphite, respectively. They are then 
combined by taking 1/3 of the first contribution, and 2/3 
of the second, respectively. This procedure is exact in the 
electrostatic limit for spheres, and is an excellent approx- 
imation in the Rayleigh limit (27ra;aoq < 1, where a eq is 
the radius of the equal- volume sphere) . It is also expected 
to be a good approximation for more complicated shapes 
in the Rayleigh limit provided one assumes that the c- 
axis of graphite is randomly oriented with respect to any 
symmetry axis that might be present within the particles. 

There is no consistent report of the interstellar UV 
feature being polarized (i.e. the particles are not partially 
oriented elongated particles). Therefore, in this paper, we 
assume that the particles are randomly oriented. 

It has been confirmed (Witt 1989; Witt et al. 1992; 
Calzetti et al. 1995) that the extinction in the vicinity of 
the interstellar UV feature is consistent with pure absorp- 
tion (i.e., carrier with size in the Rayleigh limit). This fact 
is exploited to simplify the orientational average by taking 
only three mutually perpendicular orientations of the inci- 
dent electric field with respect to the particles. This aver- 
aging procedure is also exact in the electrostatic limit and 
is a good approximation in the Rayleigh limit. An added 
advantage of the Rayleigh limit is that one needs not worry 
about size distributions of grains since extinction is inde- 
pendent of size (except perhaps indirectly through grains 
of different sizes possibly having different shapes). 

It is of interest to see how the width, shape, and peak 
position of the graphite UV feature change in the case of 
ensembles of agglomerated spheres to simulate clustered 
or irregularly shaped interstellar grains. Such a configura- 
tion can be handled by a multiple Mie sphere code (Mack- 
owski 1991; Rouleau 1996), provided one assumes that the 



spheres, though touching, do no interpenetrate each other, 
and are isotropic. Therefore, we have performed multiple 
Mie sphere computations to calculate extinction cross sec- 
tions per unit volume, C/V, for various simple arrange- 
ments of spheres and plausible clustering models. 

The two types of clusters considered are a cluster- 
cluster agglomerate (CCA), and a compact cluster. A 
cluster-cluster agglomerate is constructed from hierarchi- 
cally joining together clusters containing the same num- 
ber of spheres, N, in an N = 2, 4, 8, 16, . . . progression 
(see, e.g., Stognienko et al. 1995). These are rather fluffy 
and irregular in shape (in the case N = 16 considered 
here, the filling factor with respect to a sphere enclosing 
the cluster is about 0.04). Another alternative is to use 
a particle-cluster agglomerate (PCA) in which individual 
spheres are added to the cluster with a sticking probabil- 
ity of unity. For the small number of spheres used here, 
both CCA and PCA clusters yield similar results and so 
PCA clusters are not included. A compact cluster rep- 
resents the tightest arrangement of spheres possible, with 
each sphere touching many neighbours (Rouleau & Martin 
1993; Rouleau 1996). Again, the spheres are not allowed 
to interpenetrate. The filling factor of a compact cluster is 
thus much larger than either CCA or PCA clusters (about 
0.26 for an N = 16 cluster). This type of cluster could be 
the result of a low sticking probability of the grains or a 
compaction of an initially looser agglomerate after some 
disruptive event (like interstellar shocks). 

The clusters are assumed to be in the Rayleigh limit. 
The optical constants of Draine (1985) are used for the 
computations. Qualitatively similar results are obtained 
for other choices of optical constants of graphite. A mul- 
tipole expansion order up to 6 was necessary to obtain 
convergence for the C/V results with the multiple Mie 
sphere calculations. 

There are three distinct ways to deal with the 
anisotropy of graphite in the multiple Mie sphere code. 
The first way is to use the "1/3-2/3" rule, i.e. to average 
the extinction cross sections of the agglomerates with ei- 
ther eii or e±. The second way is to randomly assign en or 
e_i_ with probabilities 1/3 and 2/3 to the individual spheres 
forming the aggregates. The third way is to average the 
components of the dielectric tensor and to use the result- 
ing dielectric function in the multiple Mie sphere computa- 
tions. These different approaches should be good approxi- 
mations if we assume that there is no correlation between 
the orientations of the anisotropic spheres within the clus- 
ters and the direction of their individual c-axes. This is 
true for the regions outside of the resonance. However, it 
turns out that the results of the three methods differ from 
each other markedly in the bump region (Fig. |l|). E.g. for 
the CCA cluster the averaged dielectric function (third 
way) leads to a much broader feature (7 = 2.20 /im _1 ) 
at smaller wavenumbers (x max = 4.48 /im -1 ) compared 
to the results obtained with the "1/3-2/3" rule (results 
for other arrangements of spheres are discussed in detail 
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Fig. 1. Cross section per unit volume, C/V, of randomly ori- 
ented N — 16 CCA graphite clusters. Three different ways (see 
text) are used to treat the anisotropy of graphite: "1/3-2/3" 
approximation (long dashed line), random assignment (solid 
line), and averaged dielectric function (short dashed line). The 
vertical line marks the interstellar UV feature peak position 



Fig. 2. C/V of a Rayleigh sphere (dotted line), a bisphere 
(dash-dotted line), a chain of three spheres (dash-triple-dot 
line), an N = 16 CCA cluster (long dashed line), an N = 16 
compact cluster (solid line), and CDE (short dashed line). 
All arrangements of spheres are randomly oriented, using the 
"1/3-2/3" approximation to treat the anisotropy of graphite 



below) and the random assignment (second way) . Further- 
more, the feature obtained with the "1/3-2/3" rule shows 
a small shoulder in the long-wavelength wing. The appear- 
ance of the shoulder is explained by the rather long linear 
chains of the strong e± "oscillators". In case of the ran- 
dom assignment the chains are interspersed by the weak 
e II "oscillators" which makes the feature less sensitive to 



the detailed cluster structure (see Sect. \i.2\ ). The averaged 
dielectric function (third way) is in any case also more or 
less moderate compared to e± regardless of the used aver- 
aging procedure. The feature position and width, however, 
certainly depend on the averaging procedure (the numbers 
given above are for the Bruggeman mixing rule; Bohren 
& Huffman 1983). Which way and which averaging pro- 
cedure have to be used to get the best approximation of 
the C/V of clusters of graphite spheres may be tested 
with extensive DDA calculations in the future, but this 
subject goes beyond the goals of this paper. In any case, 
one should note that more detailed information about the 
arrangement of basic structured units (in respect to their 
individual c-axes) from laboratory data are required to get 
a more reliable basis for the calculations (see Rouzaud & 
Oberlin 1989). The remaining discussion in this section is 
restricted to the results obtained with the "1/3-2/3" rule 
for the various sphere configurations. 

Figure | shows C/V for a single sphere, a randomly 
oriented bisphere, a randomly oriented linear chain of 
three touching spheres, a CCA cluster containing 16 
spheres, and a compact cluster also containing 16 spheres. 
Other random realizations of these clusters, also contain- 
ing 16 spheres, gave very similar results. In particular 
the shoulder observed around 4^m _1 is reproduced in all 



Table 1. Computed and fitted C/V peak position ("PEAK" 
and 2 maX) respectively) and width ("FWHM" and 7, respec- 
tively) for arrangements of graphite spheres (in pm -1 ) 

cluster PEAK FWHM z max 7 y 2 



1 sphere 


4.743 


1.21 


4.762 


1.06 


67.6 


2 spheres 


4.735 


1.48 


4.738 


1.31 


74.6 


3 spheres 


4.742 


1.73 


4.731 


1.52 


225.2 


CCA 


4.716 


1.87 


4.708 


1.79 


365.4 


compact 


4.696 


1.61 


4.680 


1.48 


34.9 


CDE 


4.268 


2.46 


4.332 


2.40 


21.4 



CCA cluster realizations. The appearance of the UV fea- 
ture looks very similar if only C(ej_)/V is plotted. Note 
that in going from a single sphere to chains of increasing 
N, the peak decreases in amplitude, while the extinction 
rises preferentially in the long- wavelength wing of the fea- 
ture. The feature of a CCA cluster is similar to that of 
a chain of three spheres, except that this extinction en- 
hancement has developed into a shoulder. Also shown for 
comparison (short dashed line) is C/V for a continuous 
distribution of ellipsoids (CDE; Bohren & Huffman 1983), 
a widely used approximate clustering model. 

Table lists the peak position ("PEAK"), full- width 
at half- maximum ("FWHM") of the computed C/V, as 
well as x max and 7 derived from a fit similar in form to 
Eq. (1), with fitting coefficients ci, 02, c 3 and C4, using 
a slightly modified version of the routine given in Fitz- 
patrick & Massa (1990). To indicate the degree of quality 
of the fit, the x 2 value (in units of 10 6 ) is also given, 
which is defined in the usual way with the weight for 
each C/V data point (sampled at uniform intervals in 
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wavenumber) set to the (arbitrary) value 1 except for the 
bump region x = 3.3-5.9 /im" 1 where the weight is set 
to 2. Apart from the CDE case (which turns out to be a 
bad clustering model for the bump), the parameter cc max 
is quite stable (^ 1.7% relative variation), whereas 7 is 
substantially increased (by 40 — 70%) in going from a sin- 
gle sphere to clusters. These features satisfy qualitatively 
the observational constraints, except for the fact that the 
peak position is wrong (~ 4.68-4.76 /im" 1 compared to 
4.60 ± 0.04 /im" 1 ). Note that there is a significant differ- 
ence between "PEAK" and x max , and between "FWHM" 
and 7. This is due to the fact that a Drude profile fit 
is applied to a profile that is not purely Drude-like, thus 
introducing a spurious linear background (coefficients ci, 
C2) which modifies both the peak position and the width 
of the feature. The fit in the case of the compact cluster 
is reasonable, which is indicated by the small \ 2 value in 
Table |l|, whereas the fit is somewhat poorer in the case 
of the CCA cluster (due to the presence of the shoulder), 
yielding a larger x 2 . 

The peak position is at a larger wavenumber than that 
observed for the interstellar UV feature. This is why one 
has to resort to either shape or size effects to shift the fea- 
ture at the position of the interstellar UV feature (see, e.g., 
Draine 1988; Mathis 1994). Discounting the disagreement 
in peak position and concentrating only on the qualitative 
aspect of the feature, one can surmise that some irregu- 
larly shaped or clustered particles can in fact broaden the 
feature without appreciably shifting its peak position, as 
required for the purported interstellar UV feature carrier. 
However, "fluffy" or elongated structures tend to enhance 
the long wavelength wing of the feature, even introducing 
additional structure in the feature (like the CCA cluster 
does). Assuming this is a general feature of CCA clus- 
ters (see Sect. ||), it then appears that fluffy particles are 
severely constrained as a possible "topology" of the inter- 
stellar UV feature carrier. A compact cluster is an appeal- 
ing alternative, since its overall spherical shape and com- 
pactness suppresses the appearance of additional structure 
in the feature, but still allows the feature to be broadened 
substantially compared to isolated spheres. The broaden- 
ing, however, appears to be insufficient to explain all the 
range of FWHM observed for the interstellar UV feature. 
A further contributor to the broadening could be an in- 
trinsic variability in chemical composition. This interpre- 
tation is consistent with the fact that the narrowest in- 
terstellar UV features appear in all sorts of environments, 
diffuse and dense alike, with an accompanying variation in 
peak position over the whole observed range. We expect a 
clustering of grains to occur primarily in denser regions. 
The ensuing broadening, however, appears to contradict 
the observation that the broadest interstellar UV features 
(observed in dense regions and thus presumably arising 
from clustering) are accompanied by a shift of a; max to 
larger values, whereas here, it shifts to smaller values (see 



Table [j]). But this does not necessarily need to be the case, 
as shown in the following section. 

4. Single-Lorentz oscillator models 

Actually, in view of the widely assumed formation mecha- 
nisms of the carbonaceous component responsible for the 
UV bump, the highly anisotropic dielectric function of pla- 
nar graphite is unlikely to be a good model of the pur- 
ported dielectric function of the bump carrier. First, the 
exposure of the carbonaceous interstellar grains to UV 
radiation, though conducive to a "graphitization" of the 
material through dehydrogenation, is unlikely to produce 
perfect graphite sheets as an end-product (Sorrell 1990). 
There must still be considerable defects in the structure, 
the topology of the grains being a collection of randomly 
oriented graphitic crystallites in a sp 3 bonding matrix 
(possibly containing some hydrogen). This structure is 
more characteristic of amorphous carbon (especially UV 
processed or annealed hydrogenated amorphous carbon 
- HAC; Fink et al. 1984; Mennella et al. 1995a,b). Fur- 
thermore, modification of the optical constants of graphite 
have already been shown to be necessary in order to repro- 
duce even the most basic properties of the interstellar UV 
feature (Mathis 1994). An intrinsic variability in chemical 
composition must also be considered in order to satisfy 
the observational constraints. 

Virtually all UV bumps observed to date can be re- 
produced extremely well by a "Drude" profile. This is 
the profile generated by a sphere whose optical proper- 
ties are characterized by a single-Lorentz oscillator model. 
Thus, in this section, we consider a single-Lorentz oscilla- 
tor model to approximate the dielectric function hypoth- 
esized for the interstellar UV feature carrier, and see how 
shape and clustering can affect the peak position, width, 
and shape of the UV bump. We assume that variability in 
the width is attributable to shape and clustering effects, 
as well as to intrinsic variations in chemical composition 
along different lines of sight. Variability in the peak po- 
sition, which is observed to be uncorrelated with width, 
is attributed to variations in mean chemical composition 
alone. 

The dielectric function of a Lorentz oscillator is given 
by e = 1 + x 2 p j\x\ — x 2 - i'yox], where x p , x , and 70, 
are the plasma frequency, the peak position, and the 
damping constant, respectively, all in units of the inverse 
wavelength (/im" 1 ). A Drude model is characterized by 
xq = and usually describes metals (or semi-metals, like 
graphite). Through these parameters, the model can sim- 
ulate the variability of the chemical composition of the 
grains and provide some physical insight into the mecha- 
nisms involved. 

Interstellar grains cannot be expected to be perfectly 
smooth spheres in all environments. For example, one may 
expect aggregation of the primary grains, especially in 
denser regions. Furthermore, the grains could be char- 
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acterized by surface roughness and/or by porosity, as 
well as chemical inhomogeneities. An interesting question 
is whether such shape and clustering effects (neglecting 
chemical inhomogeneities within a given grain) can con- 
serve the Drude-like profile that is observed, along with 
the other observational constraints. 

4--1- Combining a Lorentz oscillator model with shape and 
clustering 

The effect of shape and clustering on the scattering prop- 
erties of chemically homogeneous grains in the Rayleigh 
limit can be modelled via a spectral density, g{L), of the 
geometric factor L, where < L < 1 (Bohren & Huffman 
1983; Fuchs 1987; Rouleau & Martin 1991). The require- 
ments are that the zeroth and first moments of g{L) are 
unity and 1/3, respectively. This approach is closely re- 
lated to the "Bergman representation" which is based on 
effective medium theory in the case of a binary mixture 
(Stognienko et al. 1995). The two approaches are equiva- 
lent if one assumes that one of the components is vacuum. 
However, the spectral density g(L) can be directly inter- 
preted in terms of shape and clustering of small grains, 
whereas the Bergman representation deals with bulk ma- 
terial, and so gives only indirect information about the 
scattering properties of small grains. 

Combining a single-Lorentz oscillator model of dielec- 
tric function e = 1 + x 2 /[xq — x 2 — ijo x ] with a model 
of shape and clustering using g(L), the extinction cross 
section per unit volume can be written as 



— = 2-7T x Im 



2?r 



g(L) dL 

(e-l)-i+i 



o [x - {xl + x 2 p L)/x] + 7 g 



(2) 



where Im denotes the imaginary part. The familiar form of 
C/V for spheres is obtained from g(L) — S(L — 1/3). Note 
that C/V using this model has the form of the "bump" 
term in Eq. (1) if g(L) = S(L — L a ). In that case, the 
maximum C/V occurs close to 



[4 



■ x 2 L ] 



1/2 



(3) 



One advantage of this type of representation is that 
a statistical ensemble of grains of various shapes or clus- 
tering states can also be represented by an average g{L). 
Actually some mean g{L) and mean e could even describe 
qualitatively a statistical ensemble of grains with varying 
shape and chemical composition, as long as there are no 
correlations between the two. 

The irregular grains considered here arise from the ag- 
glomeration of spheres (i.e. the spheres can touch their 
neighbour at one point, but they are not allowed to in- 
terpenetrate). In the case of agglomerates of spheres, the 
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Fig. 3. Spectral density for N = 16 CCA clusters (dashed 
line) used in Fig. ^[ and spectral density for N = 16 compact 
clusters (solid line) used in Fig. []. The vertical line represents 
the spectral density of spheres, a delta function at L — 1/3 



spectral representation approach brings about a consid- 
erable reduction in computer burden compared to direct 
computations using a multiple Mie sphere code since the 
problem has to be solved only once. Results between the 
two approaches in the case of clusters of spheres using the 
optical constants of graphite are very similar. 

To compute the spectral density g(L) we use a code 
kindly provided by Hinsen & Felderhof (1992). This code 
computes ensemble-averaged electromagnetic interactions 
between identical spheres in the electrostatic limit. As 
stated in Stognienko et al. (1995), the spectral density 
depends on the choice of the maximum order for the mul- 
tipole expansion. We present the spectral densities com- 
puted with multipole expansion order of 9 for the N = 16 
CCA clusters and the N = 16 compact clusters in Fig. |[ 
Note that a multipole expansion order of 6 is sufficient to 
obtain convergence for the C/V results using Eq. 2 with 
the optical constants of graphite. 

For our simple Lorentz model, we assign one of the 
smallest widths derived from interstellar extinction curves, 
7 = 0.77 (for HD 93028; Fitzpatrick & Massa 1986), to 
spheres. This might not be completely correct, but the 
sphere is a useful reference shape. We assume here that 
broader widths arise exclusively from shape and cluster- 
ing effects without any variation in chemical composition. 
We choose the peak position as being the mean value, 
£max = 4.60 /im -1 . This is an arbitrary choice since the 
narrowest profiles (7 < 0.8) span almost the entire range 
of x max , from 4.58 fimT 1 to 4.62 /im -1 . Other choices of 
i max would yield qualitatively similar results, but simply 
shifted in wavenumber. These would correspond to intrin- 
sic variations of the assumed chemical composition of the 
grains. 
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A" 1 ( M rTT 1 ) 

Fig. 4. C/V, e r , and t%, for various single- Lorentz oscillator 
models and the g{L) computed from N = 16 CCA clusters. 
By order of increasing amplitude, the models are given by 
xo = 4.0, 3.6, 3.2, 2.4, and 0.0 (solid lines). The dotted lines 
are the fit using Eq. (1) 

For spheres, a pure Drude profile is obtained. Due to 
shape and clustering effects, however, the synthetic bump 
computed from g(L) may differ from a pure Drude profile. 
We wish to find out under which conditions such clus- 
tering models can produce an increase in width of the 
bump without appreciably changing its peak position or 
its Drude-like shape. 

4-2. Results for CCA clusters 

Figure ^ shows the extinction cross section per unit vol- 
ume, C/V, and the real and imaginary part of the di- 
electric function, e r and Cj, for various single-Lorentz os- 
cillator models using the g(L) computed from N = 16 
CCA clusters. The models were parameterized using xo 
and x p = [3(4. 60 2 — a; 2 ,)] 1 / 2 . This choice of x p gives rise to 



Table 2. Same as Table |l| but for single-Lorentz oscillator 
models parameterized by xo in the case of N = 16 CCA clusters 
(PEAK = :r max = 4.60 and FWHM = 7 = 0.77 for spheres). 
"MAX" (cm -1 ) is the computed peak value (C/V) ma x 



XO 


MAX 


PEAK 


FWHM 


X max 


7 


2 

X 


4.59 


2.26E4 


4.600 


0.77 


4.600 


0.77 


0.0 


4.50 


2.21E5 


4.601 


0.78 


4.600 


0.78 


0.0 


4.20 


7.64E5 


4.611 


0.92 


4.601 


0.93 


28.5 


4.00 


1.02E6 


4.629 


1.06 


4.605 


1.09 


139.8 


3.80 


1.23E6 


4.652 


1.20 


4.612 


1.27 


402.9 


3.60 


1.40E6 


4.674 


1.33 


4.624 


1.47 





3.40 


1.54E6 


4.694 


1.45 


4.639 


1.69 


— 


3.20 


1.67E6 


4.712 


1.55 


4.657 


1.91 




2.80 


1.87E6 


4.744 


1.64 


4.703 


2.35 




2.40 


2.03E6 


4.771 


1.60 


4.744 


2.65 




2.00 


2.16E6 


4.793 


1.60 


4.744 


2.52 




1.60 


2.26E6 


4.811 


1.62 


4.723 


1.78 




1.20 


2.33E6 


4.826 


1.64 


4.738 


1.44 




0.80 


2.38E6 


4.836 


1.65 


4.752 


1.31 




0.40 


2.41E6 


4.842 


1.66 


4.761 


1.25 




0.00 


2.42E6 


4.844 


1.66 


4.764 


1.24 





a peak at 4.60 /zm" 1 for spheres (L — 1/3 in Eq. ||). The 
damping constant is set to 70 = 0.77/im -1 . For spheres, 
this produces a Drude profile of width 7 = 0.77/im -1 . 

In order of increasing amplitude, results shown are for 
x Q = 4.0, 3.6, 3.2, 2.4, and 0.0 (solid lines). Also shown 
for comparison are the fits using the procedure of Fitz- 
patrick & Massa (1990; dotted lines). Note the consider- 
able amount of structure in the profiles, except for the 
weakest model shown, x$ = 4.0. The corresponding g(L) 
is shown in Fig. ^| (dashed line). A comparison between 
the two figures indicates that the C/V of models with 
the smallest xq ("strong" Lorentz models, i.e. closest to a 
Drude model) are just a convolution of g(L) with a broad- 
ened Drude profile. Thus, for "strong" Lorentz models, 
extra peaks in the profile of g(L) will translate into struc- 
ture in the C/V curves. Conversely, for "weak" Lorentz 
models (x p L <C Xg in Eq. [2] ) the profile will be close to a 
"Drude" profile, irrespective of the form of g(L). 

Table || is similar to Table |l|, but for the one-Lorentz 
oscillator models using the g(L) of CCA clusters. In the 
present case, x max and 7 are meaningless for xq 5s 3.8 
since the curves are poorly represented by a "Drude" pro- 
file (thus being inconsistent with observations of actual 
interstellar extinction curves). Again, note the significant 
discrepancies between "PEAK" and x max , and between 
"FWHM" and 7 because of the spurious fitted linear rise. 
Only x 2 values of less than 500 are listed, corresponding to 
models that are acceptable. Note that, though the curves 
are not strictly Drude-like, "PEAK" shifts to larger val- 
ues compared to that of spheres (4.60 /im -1 ) for broader 
profiles, precisely what is observed in denser interstellar 
environments. This is related to the fact that the g(L) 
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Fig. 5. Same as Fig. ^, but for a g(L) computed from N — 16 
compact clusters (solid line in Fig. |3J) 

profile peaks at a value of L larger than 1 /3 (the value of 
L for spheres, indicated by a vertical line in Fig. ||). 

The ratio of carbon locked up in interstellar grains rel- 
ative to hydrogen, Nq/Nr, is thought be be around ~ 100 
part-per-million (ppm) by number (Cardelli et al. 1996). 
Assuming the bump grains have a density of 2 g cm~ 3 , 
then only the "weakest" Lorentz models (xq <; 4.5) require 
more carbon than available (cf. the values of (C/V) max ) 
in Table ||). As expected, putting all the carbon in a sin- 
gle extinction feature requires only modest amounts of the 
element. 

Thus cluster-cluster agglomeration (i.e. fluffy interstel- 
lar grains) can reproduce the observational constraints re- 
lating to the interstellar UV feature, provided the bump 
grains are described by a single-Lorentz oscillator model 
with xq ~ 3.8-4.5 /zm . 



Table 3. Same as Table |2j, but in the case of N = 16 compact 
clusters 



Xq 


MAX 


PEAK 




X max 


7 


X 


4.59 


2.25E4 


4.600 


0.77 


4.600 


0.77 


0.0 


4.50 


2.23E5 


4.600 


0.78 


4.600 


0.78 


0.0 


4.20 


7.86E5 


4.597 


0.88 


4.598 


0.88 


12.0 


4.00 


1.07E6 


4.596 


0.99 


4.595 


0.99 


54.6 


3.80 


1.29E6 


4.600 


1.10 


4.592 


1.12 


125.7 


3.60 


1.47E6 


4.607 


1.20 


4.590 


1.25 


213.9 


3.40 


1.61E6 


4.616 


1.31 


4.590 


1.39 


318.1 


3.20 


1.73E6 


4.626 


1.41 


4.591 


1.53 


448.2 


2.80 


1.91E6 


4.650 


1.60 


4.597 


1.83 


— 


2.40 


2.05E6 


4.673 


1.76 


4.609 


2.12 




2.00 


2.15E6 


4.693 


1.89 


4.627 


2.40 




1.60 


2.22E6 


4.712 


1.99 


4.649 


2.67 




1.20 


2.27E6 


4.726 


2.07 


4.671 


2.88 




0.80 


2.31E6 


4.736 


2.13 


4.691 


3.05 




0.40 


2.33E6 


4.742 


2.16 


4.704 


3.15 




0.00 


2.34E6 


4.744 


2.17 


4.709 


3.19 





4-3. Results for compact clusters 

We also computed C/V using Eq. (2) for a g(L) corre- 
sponding to an ensemble of compact clusters containing 
N = 16 spheres. Figure || is similar to Fig. |[ but for an 
ensemble of A" = 16 compact cluster models. The corre- 
sponding g(L) is shown in Fig. || (solid line). Note that 
the g(L) of compact clusters also peaks at L > 1/3, but 
that it is more "compact" and contains less structure than 
the g(L) of CCA clusters. This translates into C/V pro- 
files that are closer to a Drude profile, except again for 
"stronger" Lorentz models (xrj SS 3.6). This emphasizes 
how tightly constrained the shape and clustering of grains 
described by the stronger single-Lorentz oscillator mod- 
els are in terms of their allowable g(L) profile. Table 
lists the same parameters as Table for the various com- 
pact clusters using these Lorentz models. Models with 
3.2 Si xq Si 4.5 satisfy most of the observational con- 
straints. Only models with xq i; 4.5 violate the cosmic car- 
bon abundance constraint. Note again the trend of larger 
FWHM's usually leading to larger values of "PEAK" , as 
required observationally. 

5. Two-Lorentz oscillator models 

The simple approach taken in the previous section - to 
model the optical properties of the interstellar UV fea- 
ture carrier with a single Lorentz oscillator - may appear 
somewhat unrealistic. It assumes that the bump carrier 
does not contribute at all to the FUV rise. Actually, along 
some lines of sight (many passing through Orion), the 
FUV rise is observed to be very small or virtually absent, 
even though the interstellar UV feature is still present, 
albeit weaker than usual. But all carbonaceous materials 
(including diamond) are expected to give rise to a peak at 
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roughly 10-15 /im -1 due to the a-a* electronic transition 
(Mennella et al. 1995a; Fink et al. 1984). Moreover, any 
plausible material in the small particle limit will produce 
strong absorption at x > 6 /mi -1 . 

To our knowledge, no physical mechanism can suppress 
the a-a* transition while keeping intact the 7r-7r* transi- 
tion responsible for the UV feature. Thus, the assumption 
that the UV bump carrier is carbonaceous in nature ap- 
pears at first glance to be difficult to reconcile with the 
very small FUV rise observed in a number of cases. In 
this section, we ignore this difficulty and assume that the 
dielectric function of the interstellar UV feature carrier is 
the sum of two Lorentz oscillators - one producing the UV 
bump, and the other producing the FUV curvature- to see 
how the conclusions of the preceding section are modified. 

One important distinction must be made here. If two 
grain populations are unrelated, then one can add their 
cross sections, like is implicit in the decomposition of 
Eq.(l). For example, adding an extra linear background 
or an extra FUV curvature will affect only their corre- 
sponding fitting parameters in Eq. (1), but will leave the 
other fitting parameters unchanged. However, if the UV 
bump and the FUV curvature arise from the same par- 
ent material, then their respective contributions to the 
total dielectric function must be added instead, and then 
only can the cross section be computed assuming a cer- 
tain grain shape. These various contributions to the di- 
electric function do not necessarily add up linearly when 
translated into cross sections, unless Ci is small. There- 
fore, adding an extra Lorentz oscillator at larger x (in e 
space) to an initial one that produces a UV bump well 
characterized by a "Drude" profile (for example, assum- 
ing a spherical shape) might actually result into a bump 
that is no longer Drude-like. As a result of this process, all 
fitting parameters may be affected in a non-trivial way. 

To investigate these effects we approximate such a di- 
electric function as the sum of two Lorentz oscillators, giv- 
ing rise to features at about 4.6 /im" 1 (the UV bump) and 
10-15 /im" 1 (the FUV bump), respectively. This guaran- 
tees that the total dielectric function of the material sat- 
isfies Kramer-Kronig relations (Bohren & Huffman 1983). 
Of course, more complicated dielectric functions are also 
possible, but then the number of free parameters becomes 
rapidly excessive, and the physical interpretation becomes 
less straightforward. 

The effects of this extra Lorentz oscillator on the UV 
bump are illustrated in Fig. ||. The UV oscillator is given 
by the parameters xqi = 4.0, x p i — 3.934, and 701 = 0.77, 
which on its own yields a C/V peak at 4.60 /im _1 as- 
suming a sphere (dotted line). The FUV oscillator is 
given by xq2 = 10.0, x P 2 = 6.833, and various widths, 
which on its own yields a C/V peak at 10.75 /mi -1 as- 
suming a sphere. This could be viewed as a very crude 
model of the a-a* peak observed for HAC in labora- 
tory extinction experiments (Mennella et al. 1995a). The 
solid, dashed, dash-dotted, and dash-triple-dot lines are 



1.4x10 




6 8 10 

A' 1 (^m H ) 

Fig. 6. Same format as Figs. ^ and ^[ but for various Lorentz 
oscillator models and for spheres. Results shown are for a sin- 
gle-UV-Lorentz oscillator model with xoi = 4.0, x v \ = 3.934, 

701 = 0.77 (dotted lines), and for models involving an extra 
FUV Lorentz oscillator with £02 = 10.0, x P 2 = 6.833, and 

702 = 0.3 (solid lines), 3.0 (dashed lines), 6.0 (dash-dotted 
lines), and 12.0 (dash-triple-dot lines) 



for 702 = 0.3, 3.0, 6.0, and 12.0, respectively (also as- 
suming spheres). The main effect of the FUV oscillator 
is to decrease the strength of the UV bump, more or less 
independently of the width of the FUV bump (keeping 
the FUV oscillator strength constant). Another effect is 
to shift the peak position of the UV bump to smaller 
wavenumbers (4.506, 4.519, 4.536, and 4.573 /im -1 , for 
702 = 0.3, 3.0, 6.0, and 12.0, respectively). An extra 
broadening also occurs, but only on the larger x wing of 
the bump. Changing the width of the FUV oscillator with- 
out changing its strength affects the FUV curvature dra- 
matically. Actually, this suggests an explanation for the 
quasi-absence of FUV curvature observed along some lines 
of sight: if the FUV bump is extremely narrow (or shifted 
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Fig. 7. C/V of Mie spheres for various radii a using the dielec- 
tric function with 702 = 3.0 in Fig. ^ (dashed lines). In order 
of decreasing strength, C/V for radius a in the Rayleigh limit, 
and a = 0.01, 0.02, 0.03, 0.04, 0.05, and 0.1 /im, respectively 

to large enough wavenumbers), then the FUV rise may be 
virtually absent (solid line) . This is also obtained for very 
wide FUV bumps, but then there is a "jump" in the linear 
component of the extinction across the UV bump (dash- 
triplc-dot line). However, Fitzpatrick & Massa (1988) have 
shown that the linear background of interstellar extinc- 
tion curves, once the fitted bump has been subtracted 
out, is very smooth across the bump, with no evidence 
of jumps or changes of slope. Therefore extremely broad 
FUV bumps are more or less ruled out as an explanation of 
the lack of a FUV curvature concurrent with the presence 
of a UV bump. 

The main effect of adding such a background to the 
dielectric function is to increase e r and, to a lesser ex- 
tent, a, shifting the plasmon condition (e r < 0) to smaller 
wavenumbers. Note how little the dielectric function of the 
UV oscillator is changed for smaller wavenumbers when 
changing the width of the FUV oscillator (keeping the 
strength constant). 

When a shape different from a sphere is assumed (ei- 
ther a CCA or a compact cluster), as a result of this ex- 
tra background, the peak position is shifted to smaller 
wavenumbers, contrary to what was obtained in the pre- 
vious section. A moderate increase in FWHM is also ob- 
tained. These results are similar to what was obtained in 
the case of graphite (Table [I]) . 

5.1. The effect of scattering on the FUV curvature 

Another mechanism that can flatten the FUV bump is 
scattering by larger grains. But scattering is likely to af- 
fect the UV feature as well (mainly shifting it to smaller 
wavenumbers - which is not observed in all interstellar 
cases), unless a size distribution could be such that it 



would affect only wavenumbers larger than ~ 4.6 fan . 
However, this is unlikely, as shown in Fig. [?|. Plotted (in 
order of decreasing strength) are the C/V Mie curves 
for spheres with radius a in the Rayleigh limit, and 
a = 0.01, 0.02, 0.03, 0.04, 0.05, and 0.1 fim, respec- 
tively, using the dielectric function with 702 = 3.0 in Fig. || 
(dashed lines). Note that the UV bump is strongly affected 
in peak position and shape for a <; 0.03 /im, even though 
the FUV peak is still quite large (with much FUV cur- 
vature). The size which produces a FUV bump that is 
essentially flat, a = 0.1 /im (lowest curve), also yields a 
UV bump that is virtually unrecognizable. Here, the UV 
peak is more strongly affected than the FUV peak because 
the dielectric function in the UV range is closer to a plas- 
mon condition (see Fig. ^). Therefore, a size distribution 
of grains (weighting these curves by the volume of indi- 
vidual grains) is unlikely to flatten the FUV rise without 
severely affecting the shape and peak position of the UV 
bump as well. 



5.2. Synthetic interstellar extinction curves 

A direct comparison between the fitting parameters of in- 
terstellar curves and some synthetic ones can be useful to 
assess the qualities of the Drude profile for the UV bump 
and the term fp\jy(x) for the FUV curvature when two- 
Lorentz oscillator models are considered. To make such a 
comparison, a few simplifying assumptions are required. 
The total optical depth is assumed to arise from two con- 
tributions, one from grains producing most of the linear 
rise between 3 and 8.5 /an" 1 (labelled I) and one from the 
grains producing the bump, the rest of the linear rise, and 
the FUV curvature (labelled b). The component I is as- 
signed to silicate grains. The lack of correlation between 
the linear rise parameters and the bump parameters sug- 
gests that the grains responsible for the UV bump do not 
contribute appreciably to the linear rise (Jenniskens & 
Greenberg 1993). Thus assigning most of the linear rise to 
silicate grains is reasonable. The different slopes observed 
in the interstellar extinction curves (see, e.g., Fig. ||) could 
be assumed to arise from different size distributions of sili- 
cate grains (with larger mean grain sizes producing flatter 
curves). For the sake of simplicity, we shall not be con- 
cerned here with the detailed modelling of this component, 
but will simply assume that its C/V can be described as 
a simple linear rise in the range 3-8.5 /im -1 . For the car- 
bonaceous component 6, a key normalization parameter is 
the carbon to hydrogen ratio (by number), Nq/Nh. 

Assuming Ay = 1.7 x 10~ 22 V H i?v (Mathis 1994) and 
using A(x) = 1.086(NiCi + N b C b ), where Nj is the column 
density of dust and Cj is the cross section of component 
j (where j = I, b), we have from the definition of k(x), 

A(x) 

k(x) — Ry—; -Rv 

Ay 
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Fig. 8. A(x)/Ay, e r , and ei, for various Lorentz oscillator 
models assuming spheres. The parameters of the models, listed 
in Table ^, were chosen to fit actual interstellar curves (solid 
lines). Models are for HD 204827 (dotted lines), HD 229196 
(dashed lines), and HD 37023 {9 1 Ori D; dash-dotted lines: 
thick — Nc/Nh = 90 ppm, thin — 45 ppm) 



1.086 



1.7 x icr 22 iv H 



C b (x) 



Rv, 



(4) 



where Ry = Ay/E(B-V) is the total-to-selective extinc- 
tion ratio in the V band. We have Nb ~ 12mnNc / (pbVb) , 
where mjj is the mass of the hydrogen atom, and pb and 
Vb are the density and the mean volume of the grains, re- 
spectively. The density of bump grains is again assumed 
to be 2gcm~ 3 . The computed cross section per unit vol- 
ume of the bump grains using two Lorentz oscillators (i.e. 
seven free parameters: Xoj, x p j, joj, j — 1, 2, and Nc/Nh) 
is fitted by 



C b (. 



x) 



= Cbl + C b 2X - 



'(x- 



2 Cb3 , „ 2 + cm/fuv (x) (5) 

x vciaxl x ) > 7" 



between 3 and 8.5 /xm . Assuming that Ci(x) is linear 
with respect to x, Eq. (4) can be rewritten as 



k(x) = RyC 



en + ci 2 x + 



C b (x) 
V h 



Ry, 



(6) 



where en = (oi + Ry)/(RyC ) - c b \, c n = a 2 /{RyC ) - 
c b2 , and RyC Q ~ 0.064iV c /iVH. 

Interstellar extinction curves along specific lines of 
sight were modelled using this two-Lorentz oscillator 
model assuming Rayleigh spheres over the whole range. 
Using the same optical constants and the same normaliza- 
tion, results for CCA clusters and compact clusters were 
also computed to see how the bump parameters changed. 
Uncertainties related to the extrapolation beyond the ob- 
served range (3-8.5 /im" 1 ) do not warrant taking the ex- 
tra complication of scattering into account. The value of 
Nc/Nh was set more or less arbitrarily to 90 ppm (i.e. 
slightly below the current derived limit) and to 45 for the 
line of sight where a weak UV bump is observed (case of 
HD 37023). Smaller values than 90 ppm for the other two 
lines of sight would require more extreme Lorentz oscil- 
lator parameters and extinction curves may exhibit some 
substructure (see Fig. and so it might be more difficult 
to obtain a good fit. 

Figure U shows A(x)/Ay for spheres using e r and 
(shown in the middle and bottom panels, respectively). 
The dielectric function was computed from the parameters 
of two-Lorentz oscillator models. These parameters, listed 
in Table ^, were chosen to reproduce (assuming spheres) 
the extinction along three lines of sight spanning a wide 
range in FUV curvature (Fitzpatrick & Massa 1988; 1990). 
The lines of sight are in the direction of HD 204827 (top 
solid line, modelled by the dotted lines), HD 229196 (mid- 
dle solid line, modelled by the dashed lines), and HD 37023 
(8 1 Ori D; lower solid line, modelled by the dash-dotted 
lines). The values of Ry were taken from Cardelli et al. 
(1989). In Table |5[ the resulting fitting parameters ob- 
tained for spheres, CCA clusters and compact clusters are 
compared to those of these interstellar extinction curves. 
As can be shown by the x 2 values (relative to C/V, not 
A/Ay, in units of 10 6 ), the fit is excellent. Surprisingly, 
the fit can actually improve in going from spheres to clus- 
ters (case of HD 204827). This is mainly due to a better 
fit of the FUV curvature, not the UV bump. In case of 
HD 37023 the fit becomes slightly worse (due to a slight 
FUV curvature mismatch) in reducing Nc /AT H from 90 to 
45 ppm. Both oscillators are shifted to smaller wavenum- 
bers in case of Nq/Nh = 45 ppm, and, as expected, the 
fitted UV oscillator appears to be much stronger (larger 
x p i) than for the case of 90 ppm. Interestingly, 701 is even 
more reduced in case of 45 ppm. At any rate, the x 2 val- 
ues are usually smaller than in the single-Lorentz oscil- 
lator case (see entries xo = 4.0 in Tables || and ||). Thus 
it appears that an extra background under the UV bump 
(due to the FUV oscillator contribution to the UV oscil- 
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Table 4. Parameters of two-Lorentz oscillator models to reproduce interstellar extinction curves along selected lines of sight 





X l x pl 


7oi 


%02 


Xp2 


702 QiXlO 6 


Q 2 X10 6 


N C /N H (ppm) 




HD 204827 


4.100 4.200 


1.060 


8.200 


7.400 


2.900 


0.369 


0.151 




Qfl 

■JU 


2.60 


HD 229196 


4.105 3.890 


0.980 


9.700 


7.400 


2.350 


0.570 


0.105 




90 


3.12 


HD 37023 


4.407 2.425 


0.890 


10.200 


7.400 


2.350 


1.265 


-0.028 




90 


5.23 


HD 37023 


4.203 3.553 


0.868 


9.400 


7.400 


2.350 


2.522 


-0.052 




45 


D.ZO 


Table 5. Fitting parameters of synthetic interstellar extinction curves compared to actual 


ones 






Nc/Nn (ppm) 


x 




a I 


a •> 


as 


a 4 


y 2 

A. 






III./ ZUIOi t 






1 070 




1 91 Q 


3 901 


u.oyy 








sphere 


90 


4.623 


1.075 


-1.521 


1.219 


3.200 


0.882 


46.9 






CCA 


90 


4.553 


1.166 


-1.575 


1.235 


3.701 


0.871 
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lator) makes the bump profile less sensitive to shape and 
clustering. 

Unfortunately, given the already large number of free 
parameters (7), the models are not unique. But they give 
a qualitative idea of what the optical constants of actual 
bump grains might look like in various environments. Note 
that for curves exhibiting less FUV curvature, the second 
Lorentz oscillator must be shifted to larger x, otherwise 
the fit can be very poor in the FUV portion of the ob- 
served range. The FUV curvature observed in the direc- 
tion of HD 204827 is so large (a 4 = 0.899) that no FUV 
Lorentz oscillator with xq2 — 10/im _1 can reproduce it 
(these tend to give 04 ^ 2.2). These large curvatures can 
only be generated by bringing the FUV oscillator closer to 
the UV oscillator (xq2 ~ 8 /im _1 ). The fit is slightly poorer 
in that case for spheres since the tail of the FUV Lorentz 
oscillator significantly differs from the shape of this cur- 
vature. The bump peak position, x max , is sensitive mainly 
to Xq\ (a 1% variation in ceqi produces a 0.8% variation in 
#max ~~ which is a substantial part of the observed range of 
variation of x max ), and to a lesser extent, x v \. The bump 
strength 0,3 is sensitive to xq2 and x p i, whereas the cur- 
vature 04 is sensitive to X02 ■ Table || also shows that clus- 
tering, though still broadening the UV bump compared 
to spheres, systematically shifts the peak position of the 



bump to smaller wavenumbers, which is contrary to what 
was found for single-Lorentz oscillator models (and con- 
trary to observations). Thus, the inclusion of a background 
dielectric function (arising from the oscillator associated 
with the FUV peak) inhibits the conditions under which 
a shift to larger wavenumbers is possible. Therefore, if the 
optical constants of Fig. || are representative of the in- 
terstellar UV feature carrier, then intrinsic variations in 
chemical composition, rather than clustering, may be re- 
sponsible for this shift. 

6. Physical mechanisms of relevance 

In this section we discuss some of the laboratory studies 
examining the UV properties of some carbonaceous mate- 
rials of astrophysical relevance and assess their usefulness 
in interpreting the results of the preceding sections. 

Fink et al. (1984) have derived q of HAC thin films for 
various annealing temperatures, T a . As T a was increased, 
the 7T-7T* peak (associated with sp 2 bonding) shifted to 
smaller wavenumbers and its strength increased. Equat- 
ing a larger T a with less hydrogenation, this means that 
we expect xqi to shift to smaller values and x p i to in- 
crease as the hydrogen content decreases. Figures || and || 
show qualitatively that trend. On the other hand, as T a 
was increased, the a-a* peak (observed for virtually all 
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carbonaceous materials, including diamond) initially at 
~ 10/itm -1 , shifted to larger wavenumbers (~ 11/Ltm -1 
at T = 540°C), but its strength stayed roughly the 
same. Thus from this, one expects qualitatively less hy- 
drogenated material (larger x p i) to exhibit larger xq2 (and 
thus, smaller curvature 0,4) . This appears to be contra- 
dicted by the results of Fig. ||, though the comparison may 
not be relevant since the widths observed in Fink et al. 
(1984) are much broader than those considered here. 

A dehydrogenation study of small HAC grains in ex- 
tinction has been carried out by Mennella et al. (1995a). 
The HAC grains "as produced" show no UV bump, just 
a UV rise which is the small x tail of the a-a* electronic 
transition of carbon (observed at about 10.75 /1m" 1 ). As 
the annealing temperature is increased from 250° C to 
800°C, a weak UV bump develops at around 5 /urn -1 and 
gradually shifts to smaller x, gaining in strength. This is 
compatible with the above trend observed in e;. When 
all hydrogen is lost (at around 800° C), the very broad 
UV peak centered around 4^m _1 resembles that of arc- 
evaporated soot produced in an inert gas atmosphere (and 
UV features associated with hydrogen-poor circumstcl- 
lar environments; Blanco et al. 1995). These features are 
much weaker and much broader than the interstellar UV 
bump. 

Simple dehydrogenation, by itself, cannot be the only 
process responsible for the interstellar UV feature since it 
violates the observed stability of the bump peak position. 
This stability endures despite wide variations in tempera- 
ture, density, and UV flux arising from various interstellar 
environments. 

Very recently, Mennella et al. (1996) have reported a 
UV feature in small carbonaceous grains falling close (at 
2150 A) to the peak position of the interstellar UV bump. 
This feature was produced by subjecting small HAC grains 
to UV radiation (corresponding to doses about 7 times 
less than typical in the diffuse interstellar medium) . More 
importantly, this feature was shown to be stable in peak 
position when subjected to various doses of UV radiation. 
However, the precise mechanism causing this observed sta- 
bility is still sketchy at the moment. After the UV process- 
ing, the grains still contained a considerable amount of 
hydrogen (about 0.3 relative to carbon by number, about 
half the value found in the starting material). Unfortu- 
nately, the UV features produced are considerably broader 
(<; 4 /jm" 1 ) than the interstellar UV feature. This could 
be due to the extreme clustering observed for the grains 
and possibly size effects (since the individual grains have 
a radius of about 50 A) . This laboratory cosmic dust ana- 
logue nevertheless looks extremely promising in providing 
a better understanding of the optical properties of the in- 
terstellar UV feature carrier. 



7. Conclusion 

In this paper, constraints on the possible clustering state 
and dielectric function of the interstellar UV feature car- 
rier were derived using various models. 

Compact clusters of graphite spheres were shown to 
satisfy qualitatively the observational constraint that the 
UV bump peak position is very stable while its width ex- 
hibits a wide range of values. The main drawbacks of this 
model were that the peak position was wrong, and that it 
did not allow for variations in chemical composition. The 
latter is essential to account for the fact that variations 
in peak position and width are not correlated, except for 
the widest bumps (observed in dense interstellar regions 
where clustering of the grains is expected). 

To remedy this, a single-Lorentz oscillator model in 
conjunction with clustering was considered. It was shown 
to satisfy the observational constraints, including the (cor- 
rect) correlated shift in peak position observed for very 
wide bumps. 

But since the UV bump is usually assigned to some car- 
bonaceous material, then the carrier of this feature must 
also produce a FUV feature, the tail of which would be 
the FUV curvature observed along most lines of sight. A 
relatively simple two-Lorentz oscillator model was shown 
to reproduce the extinction along various lines of sight for 
a wide range of FUV curvatures. However, in this model, 
clustering shifted the peak position to smaller wavenum- 
bers for the widest bumps, contrary to what is observed. 
This model allowed a derivation of some optical constants 
of bump grains along specific lines of sight, although the 
solution was not unique. Thus, more measurements of the 
FUV part of the spectrum of carbonaceous interstellar 
dust analogues would be most desirable to further con- 
strain the model. 

Acknowledgements. The computations in Sect. ^ were per- 
formed on the CRAY T-90 of the Zentralinstitut fur Ange- 
wandte Mathematik, in Julich, Germany. We acknowledge 
the constructive comments of the referee of this paper, B.T. 
Draine. 

References 

Blanco A., Fonti S., Orofino V. 1995, ApJ, 448, 339 

Bohren C.F., Huffman D.R. 1983, Absorption and Scattering 

of Light by Small Particles, John Wiley, New York 
Calzetti D., Bohlin R.C., Gordon K.D., Witt A.N., Bianchi L. 

1995, ApJ, 446, L97 
Cardelli J.A., Clayton G.C. 1991, AJ, 101, 1021 
Cardelli J.A., Savage B.D. 1988, ApJ, 325, 864 
Cardelli J.A., Clayton G.C, Mathis J.S. 1989, ApJ, 345, 245 
Cardelli J.A., Meyer D.M., Jura M., Savage B.D. 1996, ApJ, 

467, 334 
Draine B.T. 1985, ApJS, 57, 587 
Draine B.T. 1988, ApJ, 333, 848 
Draine B.T., Lee H.M. 1984, ApJ, 285, 89 
Draine B.T., Malhotra S. 1993, ApJ, 414, 632 



14 



F. Rouleau et al.: Constraints on the properties of the 2175A interstellar feature 



Fink J., Muller-Heinzerling Th., Pfliiger J., et al. 1984, Phys. 

Rev. B, 30, 4713 
Fitzpatrick E.D., Massa D. 1986, ApJ, 307, 286 
Fitzpatrick E.D., Massa D. 1988, ApJ, 328, 734 
Fitzpatrick E.D., Massa D. 1990, ApJS, 72, 163 
Fuchs R. 1987, Phys. Rev. B, 35, 7700 
Henrard L., Lucas A.A., Lambin Ph. 1993, ApJ, 406, 92 
Hinsen K., Felderhof B.U. 1992, Phys. Rev. B, 46, 12955 
Jenniskens P., Greenberg J.M. 1993, A&A, 274, 439 
Mackowski D.W. 1991, Proc. R. Soc. Lond. A, 433, 599 
Mathis J.S. 1994, ApJ, 422, 176 

Mathis J.S., Rumpl W., Nordsieck K.H. 1977, ApJ, 217, 425 
Mennella V., Colangeli L., Blanco A., et al. 1995a, ApJ, 444, 
288 

Mennella V., Colangeli L., Bussoletti E., et al. 1995b, ApJS, 
100, 149 

Mennella V., Colangeli L., Palumbo P., Rotundi A., Schutte 

W., Bussoletti E. 1996, ApJ, 464, L191 
Rouleau F., 1996, A&A, 310, 686 
Rouleau F., Martin P.G. 1991, ApJ, 377, 526 
Rouleau F., Martin P.G. 1993, ApJ, 416, 707 
Rouzaud J.N., Oberlin A. 1989, Carbon, 27, 517 
Snow T.P., Buss R.H., Giha D.P., Swings J.P. 1987, ApJ, 321, 

921 

Sorrell W.H. 1990, MNRAS, 243, 570 

Stognienko R., Henning Th., Ossenkopf V. 1995, A&A, 296, 
797 

Witt A.N. 1989, in IAU Symp. 135, Interstellar Dust, ed. L.J. 
Allamandola and A.G.G.M. Tielens (Dordrecht: Kluwer), 
87 

Witt A.N., Petersohn J.K., Bohlin R.C., et al. 1992, ApJ, 395, 
L5 



This article was processed by the author using Springer- Verlag 
LT E X A&A style file L-AA version 3. 



